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Abstract 

We consider a structural credit model for a large portfolio of credit risky assets where the 
correlation is due to a market factor. By considering the large portfolio limit of this system 
we show the existence of a density process for the asset values. This density evolves according 
to a stochastic partial differential equation and we establish existence and uniqueness for the 
solution taking values in a suitable function space. The loss function of the portfolio is then 
a function of the evolution of this density at the default boundary. We develop numerical 
methods for pricing and calibration of the model to credit indices and consider its performance 
pre and post credit crunch. Finally, we give further examples illustrating the valuation of 
exotic credit products, specifically forward starting CDOs. 

1 Introduction 

The rapid growth of the credit derivatives market from 2000-2007 led to the development of 
increasingly complex credit instruments requiring new mathematical models for pricing and risk 
management. The subsequent contraction due to the credit crunch has placed even more emphasis 
on the importance of understanding the risks involved in dealing with complex credit products. 
Our aim in this paper is to extend standard large portfolio credit models by introducing dynamics 
and working with the infinite dimensional limit. Although this model has shortcomings (as in- 
herent in the underlying structural model) , we provide a mathematical basis for the development 
of more realistic extensions. 

The two natural approaches to credit modelling that have been extensively developed are the 
structural approach and the reduced form approach, and each has been extended to the portfolio 
setting in a variety of ways. We consider a dynamic large portfolio model obtained by taking 
the large portfolio limit of a multidimensional structural model. By taking this limit we obtain a 
stochastic partial differential equation which models the evolution of the value of a large basket 
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of underlying assets. The key quantities for multiname credit are then certain functions of the 
solution of this stochastic partial differential equation. 

Our motivation for the development of our structural evolution model came originally from the 
lack of dynamics in the credit market's standard pricing methodology. This absence of dynamics 
made pricing some structured credit instruments very difficult and credit market developments 
since mid-2007 have further exposed the limitations of existing approaches. There is still a need for 
a new generation of models to enable a better fitting as well as understanding of the risks inherent 
in some of the more complex products. For instance the existence of 5, 7 and 10-year index and 
bespoke tranches requires a model that can fit the entire correlation skew term structure, not just 
the correlation skew for a given time horizon. Also the introduction of forward starting tranches, 
options on tranches and STCDOs with trigger features requires information on the dynamics of 
spreads and information on the timings of default for their pricing. By investigating the behaviour 
of our simplified model, we are able to gain an insight into which aspects of dynamic models are 
important for the pricing of more exotic structured credit products. This information can then 
be used to help guide future model development. 

Our model follows a bottom-up approach in which the individual entities in a credit basket are 
modelled. This approach (whether structural or reduced form) has been widely used, primarily as 
a result of the introduction of copulas and the subsequent conditionally independent factor (GIF) 
models. These models allow the problem of specifying the marginal distributions and the market 
co-movements to be separated and through the choice of specific copulas has led to simple, easy 
to implement and computationally efficient techniques for pricing credit products. 

However, as the portfolio credit market expanded it became clear that these models were unable 
to cope with some of the new instruments. Copula and GIF models have no dynamics to speak of; 
nowhere is it specified how their parameters or underlyings evolve. Furthermore, they only model 
expected defaults within one time period (making copula parameters time dependent leads to 
prices that are not arbitrage free). For instruments such as coUateralised debt obligations (GDO) 
this is not an issue as they are essentially one period instruments, but for those with stronger 
timing features this is not acceptable. These two points make it impossible to price dynamic 
instruments such as options on GDOs and very difficult to price multi-period instruments such as 
forward starting GDOs. Thus our purpose is to develop a relatively simple dynamic extension of 
a GIF to the large portfolio setting. 

An alternative multi-asset route is a top-down approach where the joint default distribution 
is modeled directly without regard to the single name market. The correlation is an inherent 
property of the quantity being modeled and thus does not need to be specified. Using the top- 
down approach, frameworks similar to that of the HJM interest rate models have been developed 
for the joint loss distribution. 

Although many of the exotic credit instruments have traded infrequently, especially post the 
credit crunch, their introduction highlighted the need for a more sophisticated approach to port- 
folio credit modelling. There is a large and rapidly growing literature in this area, so we only 
mention a few papers [5], [H], [13], [3], [33]. Top-down approaches include the Markov chain 
model in [40] and the models of [4], [12] and [10]. Reduced- form approaches have been extended 
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to more than one issuer via correlated stochastic parameters. A relatively tractable example is the 
intensity-gamma model by ; another is the afRne jump diffusion model of [33] . [33] provides 
an overview of some of the main bottom-up approaches. 

1.1 Structural models 

Our model falls into the class of multi-dimensional structural models and we take the approach 
of modelling the empirical measure of the asset prices in the basket when the underlyings have 
dynamics linked through a factor model. The pricing of CDOs is then a function of the limit of 
the empirical measure of the large basket. 

Structural models are based on the premise that when a company's asset value falls below a 
certain threshold barrier a default is triggered. The first model of this type was introduced by 
[3Tj and then extended by [5]. To date, there are many variants of this model but the basic type 
is as follows. Let At be the asset value of a company whose evolution is governed by 

= fidt + a dWt , 

where /x is the mean rate of return on the assets, a is the asset volatility and Wt is a standard 
Brownian motion. If we denote the default threshold barrier by Bt we define the distance to 
default, Xt, as 

Xt^-{\ogAt~logBt). (1.1) 
a 

The event of a default by time t is now expressed as the event that X hits before time t. 

Structural models are appealing due to their intuitive economic interpretation and the link they 
provide between the equity and credit markets. They introduce spread dynamics and allow market 
participants to hedge spread risk with the underlying equity of the reference entity. Defaults are 
endogenously generated within the model and recovery rates do not need to be determined until 
after a default occurs. 

There are however downsides that affect the practical applicability of structural models. Due 
to the diffusive nature of the asset process, and the assumption of perfect information regarding 
asset values and default thresholds, any credit event generated by the model is predictable. The 
immediate consequence is short term credit spreads that are near zero: a fact contradicted by 
empirical evidence. Extensions that try to address these issues include CreditGrades"^*^ described 
in [TS] , as well as [TT] , [55] , [33] , [TO] and [B] . As structural models are extended in these ways their 
analytic complexity increases dramatically. Credit spread prices can then no longer be expressed 
in closed form and numerical methods must be employed for pricing. Another downside is that 
calibration of the model parameters is not a straightforward exercise. 

Due to the popularity enjoyed by GIF and copula models, multidimensional structural models 
have typically received less attention; as a result, the literature on this subject is relatively sparse. 
The first authors to incorporate default correlation into first passage models were [l^ and |20j . 
The former extended the Black-Gox framework to include correlated asset value processes, with 
hitting times being calculated from a time dependent barrier in closed form for two risky assets. 
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[2D] followed Zhou's approach and moved to a higher dimensional space but had to sacrifice the 
analytic results. In [21] the asset value processes for a multi-dimensional structural model are 
correlated via a set of common factors. In this setting piecewise default barriers are calibrated to 
match market prices and Monte-Carlo simulation is used to value single tranche CDOs (STCDOs). 
Other recent papers using a structural approach include [16] , [17] , [8] and [7] . We aim to develop 
a model which can allow pricing of exotic options on CDO tranches and note that there has been 
some discussion of such products in [52] , [23] . 

1.2 The SPDE model 

The starting point for our model is very similar to that used in [21] . We will develop a simple 
model in this paper in which all assets have the same constant volatility and are correlated via a 
single market factor. A more general version, in which the volatility and correlation are functions, 
can be found in [33]. Let (f2^, J^^,P^) denote a probability space for a market consisting of N 
different companies whose asset values At at time t evolve under the risk neutral measure 
according to a diffusion process given by 



dAl = rA\ dt + a^l - pA^ dWl + a^A\ dMt, i = l,...,N (1.2) 

up until the hitting time of a barrier or the horizon T. We assume and Alt are Brownian 
motions satisfying 

d[Wi,Mt] =0 Vi 



and 

d\wi,w^' 



= Sij dt, 



where we have written [., .] for the quadratic covariation and will use [.] for the quadratic variation, 
and CT > is a constant and p G [0, 1) is the constant correlation. Note the co-dependence between 
the asset processes is provided solely by the Brownian motion Mt which can be thought of as a 
market wide factor influencing all of the assets. 

Thus we can write (|1.2p in terms of the distance to default process XI = (InAJ — In _B')/cr, 
with constant barrier i?*, as 



dXl = ndt + ^/T^dWl + y^dMt, t<T^, 
XI = 0, t>T^, 



XI = x'> 0, 



(1.3) 



= M{t:Xl=0}, 



for i — 1,2, N, where fi = {r — ^a'^)/a. 

It does not matter how we label our assets so make the following assumptions. We will assume 
that {Xq, X^} is a family of exchangeable, [Cb, oo)-valued random variables with E(Xq) < oo, 
where the constant Cb > 0. We assume that this initial distribution is independent of {W^} and 
M. 
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By construction we see that our system extends to an infinite system as iV — > oo and we 
will show that there is a limit empirical measure whose density satisfies an SPDE. We will write 
(51, T ^ P) with associated expectation operator E for the limit probability space containing the 
full infinite asset value model. 

In order to state our main mathematical result we need some further notation. Let (fl*^ J^^^, P^) 
be a probability space supporting a one-dimensional Brownian motion {Mt,J^t)- Let Q^^ de- 
note the cr-algebra of predictable sets on 51*^ x (0, oo) associated with the filtration T^' and 
H\{Q, c^)) = {/ : / e L2((0, oo)), /' e ^^((o, oo))}, where ^^((o, oo)) = {/ : j;^ pdx < oo}. We 
write L'^in^'^ x{0,T),g^' , H\{0,oo))) = {/(w,i,.) : /(w,t,.) e H\{0, oo))J{uj,t, .) is J-f -measurable, 
E*^ /cT < °°}- We also write for a Dirac measure at the point x. 

Let VN^t denote the equally weighted empirical measure for the entire portfolio given by 

1 ^ 

^^:^^nT.^XI- (1-4) 

1=1 

Theorem 1.1. The limit empirical measure Dt — limAr_>.oo i^N,t exists and is a probability measure 
with two components, i>t — LfSo + Vf The measure Vt is a measure on (0, oo) with density v{t, x), 
which IS the unique solution m L'^{n'^' x (Q,T),g^\ H\{0,oo))) of the SPDE 

dv = (r - ^(T^) dt + ^v^-j; dt - y/pv^ dM{t), 

(1.5) 

v{Q,x) = vo{x), w(t,0) = 0. 
The weight of the Dirac mass at is the loss function 

/•oo 

Lt = 1 — / v{t, x)dx. 
Jo 

The price of the credit products that we consider are functions of the loss function Lf. There 
is no analytic solution for this SPDE, though it can be viewed as the Zakai equation for a filtering 
problem, and thus we require numerical techniques for its solution. One natural approach is just 
to use a Monte Carlo technique to simulate the whole basket, and for small sizes of basket this 
would be a natural approach. However, as the basket size increases, the numerical solution of the 
limit SPDE becomes more computationally efficient and we discuss this in our simplified setting. 

An outline of the paper is as follows. We begin with a description of the mechanics and basic 
valuation methods of synthetic collateralised debt obligations in Section [2] in order to provide the 
necessary background for later sections. The mathematical core of the work is in Section |3] where 
we develop our infinite dimensional model for portfolio credit starting from a multidimensional 
structural model and prove Theorem ll.il We make strong assumptions with the aim of delivering 
a relatively simple, tractable model that encapsulates the information required to calculate the 
loss distribution for a portfolio of risky assets. The aim in Section |4] is to develop a suitable 
numerical scheme for solving the SPDE. Section [5] discusses the calibration and performance of 
the model when pricing tranches of the iTraxx before and after the credit crunch. 
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2 Collateralised debt obligations 



Collateralised Debt Obligations (CDOs) are securitized interests in pools of credit risky assets. 
These assets can include mortgages, bonds, loans and credit derivatives. The CDO repackages the 
credit risk of the reference portfolio into multiple tranches that are then passed on to investors. 
Prior to the 'credit crunch' the synthetic CDO, credit indices and single name Credit Default Swap 
(CDS) market together made up the majority of the total traded notional in the credit derivative 
market. However the index tranche market is currently the only area that is still active. The 
bespoke CDO business has yet to return although there are a few signs of activity. 

Although there are many different types of CDO, here we will be focussing on what is known 
as a synthetic CDO i.e. one whose collateral pool consists entirely of credit default swaps. It 
is possible to trade single tranches within a synthetic CDO without the entire structure being 
constructed. In this case the two parties of the transaction, the protection buyer and protection 
seller, exchange payments as if the CDO had been set-up. The performance of this single tranche 
CDO (STCDO) is dependent on the number of defaults that occur in the reference portfolio 
during the lifetime of the contract. 

Each tranche is defined by two points that determine its place within the capital structure: 
the attachment point and the higher valued detachment point. These are usually expressed as a 
percentage of the total portfolio notional. The tranche notional is defined as the difference between 
the attachment and detachment points. When losses arc incurred (the loss is the notional of the 
defaulted entity corrected for recovery), and the cumulative loss in the collateral pool is between 
the attachment and detachment point, the seller pays the buyer an amount equal to the loss 
incurred within the tranche. The tranche notional is then reduced by this amount. This means 
that when the cumulative loss exceeds the detachment point the tranche notional is zero. In 
return for this protection, the buyer pays a quarterly premium based off a fixed spread and the 
outstanding tranche notional. 

Say we have N entities in our reference credit portfolio each with notional A^o- We define the 
total loss if on the portfolio as 

N 

Lt = ^La{n<t}. (2.1) 

i=l 

where Li = No{l — Ri), Ri and Tj are the recovery rate and default time of the i-th entity 
respectively. If we assume the recovery rate is the same across all credit entities and equal to a 
value R then we can write 

N 

Lt = No{l-R)J2Mn<t}- (2.2) 
The outstanding tranche notional, Zt, of a single tranche within a synthetic CDO is given by 

Zt = [d-Lt]+ -[a-Lt] + , (2.3) 

and the tranche loss as 

Yt = [Lt - a]+ - [Lt - d\+, (2.4) 
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where a is the tranche attachment point and d is the tranche detachment point. 

As for a Credit Defauh Swap (CDS) the value of a STCDO is given by the difference between 
the fee leg and the protection leg. The protection buyer pays a regular fixed spread on the 
outstanding notional of the tranche. We denote the payment dates by T^, 1 < i < n, the intervals 
by 5i — Ti — Ti_i and the value of a bank account at time t by b{t). Then the value of the fee leg 
is given by 



where the expectation is with respect to a suitable pricing measure. The protection seller only 
makes payments to the buyer when the tranche incurs losses, and the value of this payment is 
equal to the change in the tranche loss Yt. However, we can express the value of the protection 
leg in terms of the outstanding tranche notional Xt as follows 



assuming that the losses are paid at the coupon dates. As in a CDS contract the par spread s of 
the tranche is chosen to make the initial value zero hence is calculated as 



From p.5|) and (12.61) we see that the key to finding the par spread is obtaining the distribution of 
the outstanding tranche notional; from (|2.3|) . this is equivalent to finding the distribution of the 
loss Lt. As all portfolio credit derivatives are essentially options on this loss variable the heart of 
every multiname credit model is determining its distribution. 

3 An infinite dimensional structural model 

Our aim in this section is to establish Theorem ll.il We will begin by describing the system 
()1.3|) by a measure valued process and showing that there is a limit empirical measure for the 
infinite system. We then proceed to establish its behaviour near before proving that its evolution 
can be captured by an SPDE. 

3.1 The limit empirical density 

Recall the equally weighted empirical measure for the entire portfolio is given by 




(2.5) 




(2.6) 



s = 



(2.7) 




We can write this as 



7 



where 

N , AT 



i=l 1=1 



Note that Z/jv,t is a loss function in that it is the proportion of companies that have defaulted by 
time t. 

Let IR+ = [0, oo). We write V{R+) for the set of probability measures on 1R+ and P(Ck+ [0, oo)) 
for the set of probability measures on Cr_^ [0, oo) where the topology is always that of weak 
convergence. We write Cp(R^)[0, oo) for the continuous 7-* (7?.-|-)- valued functions on [0, oo). 

Theorem 3.1. There exists a Cp(R_^)[0, oo)-waZMed random variable v such that 

1 ^ 

Vt = lim i>Nt= lim ir/^Sxi, P- a.s. . 



We also have a decomposition for the limit into two subprobability measures 

i>t = Lt6o + vt. 

Proof. Let us denote the system with the same dynamics but without default by {Xl). Then 

Xi = Xil^ :=F(x:,0<s<t). 

0<B<t ^ ' \ / 

Since F is independent of i, in order to show that {X*} is exchangeable in Cr[0, oo) we only need 
to show {X'} is exchangeable in Cr[0, oo). 



Since XI = + /xf + y'l — P^l + y/P^^t for all t, and {X^} is an exchangeable family, we 
have that {X^, X^} is exchangeable at any time t. 

We prove that for any N, {X^,...,X!^} is exchangeable in Ck_^[0,oo), the continuous non- 
negative functions on [0, oo). In fact, for any Borel sets Ai, A^ € (7r^ [0, oo), we need to prove 
that for any permutation a, we have 

P {X.i G Au 1.^ gAn]^¥ {x"('' eAu-.., X"(^> e An} . 

It suffices to choose the following Ai's: for any n € N, take Ai^i, Ai^n € -^(1^+) for i = 1,. . .,N 
with a time set = to < ti <■■■< tn, and set 



8 



Therefore we have, 



n / N \ 
n / N \ 



by the exchangeabihty of the increments of {X} , ...,X^} at any time t. Hence {X^, is 
exchangeable in Cm[0, oo). As a consequence we have {X^ , X^} is exchangeable in Cr^ [0, oo) 
and for a fixed {X} , ...,X^} is exchangeable in R+. As the system p.3p is easily extended to 
an infinite particle system, by de Finetti's theorem, see for example, pQ, 



1 ^ 



I/. = lim 

Af->+oo iV 

1=1 

exists almost surely in 'P{C-r_^ [0,oo)). 

We now need to show that the {vt,t £ [0, oo)} is a continuous process in the space of probability 
measures. We define a projection mapping 

Pt : Ck[0,c5o) -^R 



by setting, for any Y. e Cr[0,oo), 

Pt{Y.) = Yf 

Then define vt ■■= v o P^^ ^V{^). We first show that 



1 ^ 

Vi = lim — > Svi . 



1=1 

To establish this we denote 

N N 

i=l 1=1 

where 9^ converges weakly to 6 and 6 exists in 'P(IR) almost surely by the exchangeability of {XI} 
at any time t. For any h E Ch(M), the collection of all the bounded and continuous functions on 
M, we have 

h{x)9{dx) — lim / h{x)6N{dx). 
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Define 

N 



1=1 

then 

and UN converges weakly to i> in P(Cr[0, oo)). Tlius 

/ h(x)0{dx) = lim / fcoP-i(a;)(ajv oP-i)(da;) 
ko p-^{x)9 o p-^{dx) 
h{x)vt{dx). 

Therefore Dt — 9 = limjv^+oo Y^^=i ^xj ^ 'P'(R)- Next we want to show that vt € C-p(R) [0, cx)). 
By definition it suffices prove that when i„ — >■ Iq, we have vt^ — ^ S'fp weakly in P(IR), i.e., we want 
to show that for any open set U G /S(K), lim inf „_j.oo vt„{U) > vtg{U) ['13 , Theorem 3.3.1]. In 
fact we have, by continuity of Y and Fatou's Lemma for sets, that 

vt,{U) =voP-\U)=D{{Y.\Yt,eU}) 

/ CX3 OO \ 

U [MnYt.eu}] 



/ oo oo 
\n=l k= 



< lim inf P({K|yf, e [/}) = liminf i>t„(?7). 

Therefore, the process {vt : i G [0,oo)} exists almost surely in Cpfg) [0, oo). 

The decomposition follows from the decomposition for N companies. We then define Lt — 
i/t({0}) and to be z/j restricted to (0, oo). □ 

For a measure C.t and integrable function (f) we write 

{q^Xt) = f ^(xXtidx). (3.1) 



Let C := {f : f € C^(0, oo), /(O) = 0,\imx-^oD f{x) = 0}. Using the empirical measure (11.41) we 
define a family of processes F^''^ for (/) € C by 

1 ^ 

Pf = z^AT,,) = - ^ ^{Xi) = (0, i^jv,*) (3.2) 
1=1 



As X* = for i > Tq*, and hence 4>{Xi) = for i > Tq*, in order to apply Ito's formula to f/^'*^ 
1 



we write F^^'''' = — (l){Xl)l[t<Ti}- Thus we have 
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N,4> r-,N,, 



1 ^ /■* / 1 

vEy^ hs<Ti} [cj,'{xi)dxi + -r{xi)d[xi] 

1 ^ 

i=i 



^'{xDfids + ci>'ixi)./T^pdw: + ^'{xi)y/pdM. + -nxDds 



i-t 1 N 



rt 1 N 



•^0 4=1 i=l 

i=l 



If we define the second order linear operator A by 
we have 

Jo Jo 

We now pass to the limit by letting N oo. 

In order to determine what happens we first focus on the idiosyncratic term in 



(3.3) 



•^0 i=l 



(3.4) 



As 0' is bounded if is a martingale and, by the independence of the Wl it has quadratic variation 



[4]*= / JT:^Y.^l-p){<t>'{Xl)fds. 
■^0 i=l 



As g C there exists a constant if^ such that l^'l < K^p. Thus 



1™ / ds<^^i, 



and hence we have for any such 



^E / (1 - ^) I '^'(^^) I' ^ Vi e [o,T]. 

i=l 



Thus the random term due to the idiosyncratic component of the asset values has become deter- 
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ministic in the infinite dimensional limit and must vanish almost surely. 

Wo also note that as (j)' , (p" are bounded and i'n,s is a probability measure, we can apply the 
dominated convergence theorem to take the limit under the integrals in the other terms in p.3p . 
We summarize in the following 

Theorem 3.2. The sequence of empirical measures vpf^t on (0,oo) satisfies for all 4> £ C , 

pN,^ _^ p4, ^ as N ^oo, a.s. 

The evolution of the limit empirical measure in the weak sense is given by 



vt) = {(p, vo) + / {A(j), V,) ds + / Vs) dM,, e C. 



(3.5) 



3.2 The boundary condition 

The behaviour of vt, the limit empirical measure on (0, oo), at the boundary zero is given in 
the following theorem: 



Theorem 3.3. We have 



lim = 0, a.s.. 

elO e 



Proof. By the definition of vt, properties of weak convergence and an application of Fatou's 
Lemma, we have 



1 ^ 



1 " r 

< liminf — V P <^ < £, inf XI > 
For t < Tq, integrating the system (jl.Sp from time to t, we have: 



(3.6) 



XI ^x'+fit+ y/1 - pWi + ^Mf 



Since we know that 



where Bt is a standard Brownian motion on the same probability space, we have 



XI < e, inf Xl>0 

0<s<t 



x' + nt + Bt<e, inf (x' + u,s + Bs)>0 

0<s<t 



=P"^ \nt + < e} - P"^ int + Bt<e, inf (ns + B^) < 



/27rt 



^^^l{z-x')-^LH/2-i\z\+x'f/2t^^ 
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12-Kt 

1 



J.TTt In V J 

(> 



/27rt Jo 



-{z-fj.t-x^y /2t 



1 - e" 



<^(l-e--^) 



<- 



1 2ex' 



(3.7) 



Assume e < jCs. Since we have > Cb, ii t < ""j^i ^ , then |z — — x'| > 0, VO < z < e 



and there exists > only depending on T such that 



If i > -^fn^, then 
— ImI ' 



J_g-(2-Mt-a:')V2t < 



1 1 
^ < 



Letting :— max < C^, i — 3- >, p.7p becomes 



21^1 



Cb ^ 3 

2|m 



0<s<t 



XI < e, inf > ^ < -^ex'sC'r := a;*CT£ 



27r 



(3.8) 



where Ct is a positive constant only depending on T. Thus by p.6p and p.8p we have 

^-""■-'"<C..CiMinflf.'). 



(3.9) 



Since {Xq, is an exchangeable family of integrable random variables, limjv^.+oo T^f=i 



exists and is finite almost surely. Let K = limAr_j.oo J2i=i Now p.9p becomes 

"^t((0,e))" 



E 

By Markov's inequality, for any A > we have 

■^t((0,e)) 



< KCre. 



KCts 



therefore, for the subsequence e = 



KCt 



Thus by the first Borel-Cantelli Lemma, as A > is arbitrary and also 



> 0, we must 
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have 



lin.sup^^i«^ = 0, a.s.. 



Now for any e > 0, there exists a n such that > < £ < and hence 



Unrsup MM) < H^3up ^» ^ hnrsup ^» f - 0, 

£4-0 £ n— 



(n+1)^ 



(n+1)" 



Since 



'^*((0,e)) 



> 0, therefore 



,(i,0):=hmM^ = 0, a.s.. 



□ 

Therefore, if there is a density for the empirical measure, it wiU satisfy a Dirichlet boundary 
condition. 

Next we give an estimate on E[(i/f ((0, e)))'^] which will be needed later. In order to do this 
we require an estimate for the distribution of the first passage times of two correlated Brownian 
motions, and the Brownian motions themselves. We use a transformation to independence and 
the formula derived in 1231. 



Lemma 3.4. Let and be two correlated Brownian motions with constant correlation \ g\ < 1, 

such that 



Bq ~ ai > 0, Bq — 02 > and law Vb- Then there exists £o — \y ^~2^ 
for all £ < Eq; 



1-^2 



Q<Bl < £, inf B] >0, 0< B^ < e, inf > I < CtE 



' 0<s<t ' 0<s<t 



-2+7 



where Ct = 2^ < 
T; and 



af+a2—2gaia2 



a ^ < 



2+5 



and Kt is a constant only depending on 



TT + tan-M -^^V^ 1 , £'>0, 



tan 



Q<0. 



(3.10) 



Therefore, if g > 0, we have ^ < a < t: and 3 < 2 + ^ < 4. 



Proof. We begin by making a transformation to obtain a two-dimensional Brownian motion with 
independent components. We follow the setup and statements in Let Bt = {Bl,B^) and 

consider the process Z = a^^B, where 



We know that Z has independent components. It is easily seen that the horizontal axis is invariant 
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under the transformation T 



defined by T(x) — a ^x, while the vertical axis is mapped 



to the line zi — — 



= Z2. 



Now the time that the first Brownian motion hits zero is transformed to the time ti which 
is the first passage time of Zt to the horizontal axis; and the time that the second Brownian 
motion hits zero is transformed to the time T2 which is the first passage time of Zt to the line 
Z2 = zitana, where < a < tt is given in (j3.10p . Moreover, in polar coordinates Zt = {Rt,&t) 
starts at the point zq given by 



I a I + 0,2 — 2gaia2 
1^ '' 



and 



TT + tan 



^0 = < 



tan 



0-1 < Qa2, 
ai = ga2, 
fli > ga2- 



It is easily verified that < Oq < a. We denote by t = min(ri,T2) the first exit time of Z from 
the wedge 

Ca = {{r cos 9,r sine) : r > 0,0 < 9 < a} C . 
If z = {rcos9, rsm9) is a point in Ca we have, by |23j . 



2r 



T >t,Zt e dz} = — e +"^0)/^* \ sm 



m:9 . wkUq 



sm ■ 



n=l 



a \ t J 



drd9. 



where ly denotes the modified Bessel function of the first kind of order v 
Using this transformation and the formula p. lip we have 

Fb \ < Bl < e, inf BI > 0, < Bf < e, inf Sf>ol 

I 0<s<t 0<s<t \ 



< 



T > t,0 < Qt < a,0 < Rt < 



1-g 



^_L,-ir^^rly.t ^ ^/„./„ f ^) drd9 

Q ta a a ' \ i I 



< 



, 00 



X! ^yiT^la I — I drd9. 



ta .In . \ t / 

By the definition of the modified Bessel function, we have 



(3.11) 



(3.12) 



r/a 



E 



1 



/TTo 



\ t J ^ m!r(m + ii2L + 1) V 2t 



m=0 

CO 



1 
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\njL]\\2t J ^ (m!)2 V 2t J 

L a J m=0 ^ ' 



oo 

y-(—) 

^ ml \2t J 



.m=0 



ft 1 



where [a;] denotes the integer part of x. Using this in p. 121) we have 



0<Bl < e, inf > 0, < < e, inf > ol 

Q<s<t 0<s<t J 



< 



< 



ta 







Ti=l L a J 



t ^ 2t ' 



If we choose Eq = '° ^ ^ , then for any e < eg we have + — Srrg > 0. Therefore we can 
find a constant Kt only depending on T such that 



-e 2t < fi^j,. 



Thus 



< < £, inf > 0, < ^2 < e, inf > 



1-Q 



1-0 



2+^ 



where Ct = 2^ q rg" Kt jzr^ is a- constant only depending on g, ai, 02 and T. 

Moreover, it is obvious that < a < tt and ^<a<7rifp>0. In the latter case we have 
3<2 + -<4. □ 

a — 

Lemma 3.5. There exists eg > only depending on p and the lower hound Cb for the {Xq}, 
such that for any rj > 0, for all e < eg we have 



where Kt is a positive constant depending on T and a is given in 
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Proof. By definition of i^t, properties of weak convergence and Fatou's Lemma 



E[(j.t((0,e)))2] <E 



^ N M 

— VVe 

N M 

—T y E 



< lini inf 

Af-i-oo,Af-i-oo NM 



{0<X,'<e,^ m^^ Xj >0, 0<Xt' <e,info<s<t >0} 



lim inf 

Af^oo.M^oo NM 



{0<Xi<e,^in{^^ Xj >0, 0<Xt' <e,info<s<t Xj >0} 



(3.13) 



Since neither of the firms i or j has defaulted by time t, we have 

XI ^x' + nt+ v/1 - pW^ + y^Mt = x'+ ^it + B]- 



XI = +iit+ yjV^pWl + VP^ft = x^ + /it + Bl, 

where B\ and are correlated Brownian motions with correlation p. 
We use the Girsanov theorem (e.g. [39]) to change the measure and set 

Z,(/i) = exp (- {B\ + Bl + Aii)) , 

which is easily seen to be a true martingale by Novikov's condition. We write P for the probability 
measure on Tt given by 

P(A) := E[UZT(Ai)]; A e Tt. (3.14) 
and E for expectation with respect to P. Thus for each fixed T e [0, cx)), the process 

{(,B\.Bl) (,B\ + Mi, + /^t), -Ft, < < < T} 

is a two-dimensional Brownian motion on (il, J^t,P), where B\ and -B^ have correlation p. 

in (IS331). We have 



We now calculate the term E 



{0<Xt'<6,info<^<t Xl>a, 0<Xt^<e,mfo<a<t Xi>0} 



= E 



< <^E 



l{0<X'<e,info<,<tX'>0, 0<Xi' <e,mfo<s<t >0} Zt{p) 

1/a 



{0<Xj'<e,info<a<t Xj>0, 0<Xi'' <£,info<e<t Xi>0} 



i/b 



where 



Ji = E 



1/a 
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J2 - 



1 



1/6 



ZT{^J), 

l/a + l/h l,a > 1,& > 1, 



by Holders inequality. 
For Ji, we have 



Ji 



< < e, inf X* > 0, < X/ < e, inf > 

0<s<t 0<s<t 



1/a 



< cc* + < e, inf + > 0, < + < e, inf a;^ + > 

0<s<t 0<s<t 



1/a 



By Lemma l3^ with g = p,ai — x\a2 — we know that there exists Eq — ^ \/^^~ ^ ^'''^ ^ 



and a as in p.2p . such that for all e < Eq we have 



' I < + < e, inf + > 0, < x' + B? «, inf + > ol < Ct£^ 



As a;* > Cr and > Cr, we have 



+ - 2pa;*xJ' > ^2(1 - p)C£ 



Thus we can choose a new Eq :— ^yj+^C^B < Cb, such that for all e < Eq we have, for all i, j, 



Ji < Cre- 



For J2 we have 



J2 - 



1 



exp 



hp 



i/b 



exp 



1 + P 



Bji -\- Bj^ — 



: exp 
- exp 



1 + P 
1 + P 



E 



exp 



bp 



i/b 



i/b 



1+P 

= Jt < 00, V6, i,j e N 



i/b 



Now we have 



Eh((0,£))2] < Ji • J2 < CtJte^, Ve < eo- 
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Now for any < r/ < ^ - 1 we can choose 1 < a = (2 + ^)/(2 + ^ - 77) < (2 + ^)/3 and hence 



We will write /J^tt/q; — 77 — l>Oso that 2 + 7r/a — 77 = 3 + /3. 

3.3 The existence and uniqueness of the density 

In order to prove our main Theorem we need to recharacterise the evolution obtain in p.Sp 
as the stochastic PDE. Thus we need the measure vt to be absolutely continuous with respect to 
the Lebesgue measure to write vtidx) — v{t,x)dx for some density v. 

We introduce some notation first. Let — L^((0, 00)) be the usual Hilbert space with L^- 
norm || • ||o and inner product (•, •)o given by ||0||o = /o°° \4>{x)\^dx and (</>, V')o = /o°° 4'{^)'>P{^)d'X. 
In the following we adapt the approach in [29] to our setting. The idea to prove the existence of 
an L^((0, oo))-density is to transform our A^((0, oo))-valued process to an iJ^-valued process, by 
convolving the measure with the absorbing heat kernel, where A^((0, 00)) denotes the set of finite 
Borel measures on (0, 00). 

For any g e A^((0,oo)) and 5 > 0, we write 



We use the same notation for the Brownian semigroup on Cf,(R+), the bounded and continuous 
functions on R+, i.e.. 



IEh((0,e))^] < Kts^+^/''-\ Ve < eo. 



where Kt is a positive constant only depending on T. 



□ 




(3.15) 



where Gs is the absorbing heat kernel in R"*" given by 





We will also need to consider the reflecting heat kernel Gg{x,y), defined by 




We write the associated semigroup as 




Then it is an easy calculation to see that 



dxGs{x,y) 



dvGs{x,y). 



(3.16) 
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It is not difficult to prove the following lemma. 

Lemma 3.6. If g e MHO, oo)) andS>0, thenTsg € H". 

We will write E if the measure vt has a density which is in . Let Zs{s) — Tgiyg, where 
1/ is an A^((0, oo))-valued solution to p.Sp . Our aim is to obtain an estimate for the i7"-norm of 
the process Zs- 

It is easy to see that Ts(t) G C for any (j) G C. Thus, replacing (f> G C hy Ts(f> in p.Sp and using 
Fubini, we have 

{Zs{t),(l))o ={Ts(j>,i^t) 

^{Tscf>,i^o)+ f {fi{Tscf,y{x) + l{Ts<l>)"{x),i^,)ds+ f {^p{Ts<t>)' {x),v,)dMs. (3.17) 
Jo ^ Jo 

The integrands can be rewritten as 

POO 

{fi{Ts<l>y{x),iys) =fi / {TscbY{x)iys{dx) 
Jo 

y ^^(^J Gs{x,y)(t){y)dy^ Vs{dx) 

=/" y (yj {dxGs{x,y))(j){y)dy^ Vs{dx) 
Applying p.l6p and Fubini we have 

(/i(T5(/))'(x),zys) ^y {-dyG'six,y))<j){y)dy^ Vsidx) 

poo / poo \ 

n G^six,y)4>'{y)dyysidx) 

pOO 

=M / (rjz.,)(y)0'(y)dy 
Jo 

poo 

^-^J■ (l){y)dy{T^Vs{y))dy 
Jo 

= - IJ'{(l^,dxTl{vs))o- 
Similarly, for the term {yfp[Ts4>)' {x),Vs) we have 

{^p{Tscj,)'{x),v,) = ~^p{cp,dxTl{v,))o. 
For the term {^{Tscf))" {x) ^ Us) we can perform the same type of calculation to see 

(i(W(x),^,) = i(0,a^T,(^,))o. 
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Therefore p.l7p becomes 



1 



(3.18) 



2 JQ 



By using Ito's formula on {Zs{s), 0)^ we have 



{Zsit),^yo={Zsio),(t>)t,+ / d{Zs{s),^yo 





\2 



^{Zs{Q),4>)i+ 2{Zs{s),4>)^d{Zs{s),4,)o+ d{{Zs{s),4>)^,{Zs{s),4,)^ 
Jo Jo 

= {Zs{0),(t>)l-2fi I {Zs{s),4>)o{<l^.d^T',{vs))ods 



Jo 

t rt 

i2r 



+ / {Zs{s),<j,)o{(b,d',Ts{i^,))ods-2^ {Zsis),<j,)o{(b,d,T^{us))odM, 
Jo Jo 

+ p /*|(0,9,TJ(^.,))opds. 







We can choose a set of e C to be a complete, orthonormal basis of and taking expecta- 
tions, we have 



^\Zs{t)\\m\Zsm\l-2^IE f {Zsis),d^mi^s))ods + E f {Zs{s),dlTs{vs))ods 

Jo Jo 



+ pE / Wd^niv^ids 

Jo 

^\\Zsm\l ~ 2fiE I {Ts{v,),d^Tl{v,))od8 + E I {Ts{u,),dlTs{iys))ods 
Jo Jo 

+ pE [ \\d^T^{i,,)\\lds. (3.19) 

We now control the integral terms on the right-hand side of p.l9p in terms of the integral of 
E||T5(z^s)||q plus some constant which goes to as (5 —> 0. 

Lemma 3.7. There exist constants C^, C2 such that for S < Sq/2 we have 

E[\-2p{Ts{,ys),d^T^{,ys))o\] < \fi\-E[\\Ts{Ws\)\\l]+C^S^ + ^e~'o/^'. (3.20) 



Proof. 

{Ts{vs),d^n{vs))o^ I Ts{vs)(x)d.,Tl{vs)(x)dx 







T5{vs){x)^j^ d.,G}{x,y)vs{dy)]dx 



Ts{vs){x)[ I {d,Gs{x,y)--^e-^^),.s{dy)]dx 



\Jo 
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nCO f-OG 

= / Ts{iys){x) / dxGs{x,y)us{dy)dx 
Jo Jo 

POO 

Jo 



°° 2 _ix±vf_x + y , , , , 

, - e 25 y^i^(ly)dx 



/ Ts{ys){x)dMi's){x)dx - / Ts{v,){x) 
Jo Jo Jo 

-| />oo |>00 

;^ / 5,[(T5(j.,)(a;))2]da;- / T^(i.,)(a;) 
^ Jo ^0 

- j^Ts{vs){x) 



2 _(H±^a; + 2/ , , 
^=e -^—Vs(dy)dx 
V2n6 o 

2 (x+y)^ X + y , , , , 
:e 25 —^iy^(^dy)dx 



'2ttS 



e ' t"' ^ ) da;. 



/27rJ 



Therefore, 

|-2/i(T5(zy,),9,TJK))o| 



< 



< 



/X / (T5(z/«)(a;))2da; + /x / / 
Jo Jo \Jo 



2 _ (£+rf^ a; + w , , , , , 
_e- .5 ^^..(rfz/)) 



2 (.x+y)^ X + y , , , , 



Now let us denote 



Pi 



Jo \Jo 



^ =e '"m' ^ ^ Vs{dy) \ dx 



/27r(5 



and investigate the bound for Pi. 



Pi = 



oo / poo 



(=^+«)2 X + y 



Vs{dy) I dx 



'2tt5 



\J0 

roo .oo .oo 2 (»+.i)2+(xW (a; + tyi)(x + ^o) 



/27r5 



Jo Jo 

oc . oo 

Vs{dyr)vs{dy2) / (^4=)2e" 

ZTTO 



^2 



Vs{dyi)vs{dy2)dx 



JO 

oo />oo 







(x+yi)^+(x+y^)-' [x + 2/i)(a; + y2) , 
^ -pi -dx 

0'' 



roo o r j_ „ I 3/i+j/2 \2 f yi-y2 \2 



/o Jo 

By changing variables using 



/27r(5 



6^ 



we have 



/•oc poo f+OO 



2 _ (v-L-V2)'^ 



^2 



Since 



1 < 



< \/2 
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when z> J Ui:tM2. ^ have 



2ttS' 



2 7 

'e ^- 



2 {yi-y2)^ 



-dz 



roo roo 



/•OO 



/o Jo Jo vZird 



(52 



2 ,2 — 



,52 



■dz 



2 2 f"" /""^ ^2 _ iyi-y2) 



\/27r(5 Jo Jo 



(52 



-1 



■ OO roo 







+ CX3 



/o V2^ 



Jo 

2 roo poo 



^{yl+yl<2z'^}Vs{dyi)Vs{dy2) 



^{y,<V2z.y2<V2z>s{dyi)Vs{dy2) 



< 



2 27 
V2{-^P'e-'^-{M{0,V2z))fdz 



Therefore, 



E[|-2Ai(r5(i^.),a,Tj(j.,))o|] 

< |M|-E[||r^(i.,)||2] + HE[Pi] 

< |/i|.E[||r,(^,)||g] + HE^y^ _e-fT_(^,((0,z)))2dz 

/•OO 2 2 ;72 

= |M|-E[||T5(z.,)||^] + |^|y^ _e-%_E[^.,((0,z))2]dz. 



By Lemma l3.5l in the last section we know that for the measure- valued solution of (I3.5p . there 
exists eo > such that for all z < eq we have 

E[iy+((0,z))2] < Ktz^+^. 

Hence we have 



E[\-2fi{Ts{iys),d,Tl{iys))o\] 

/•eo o 2 ^2 /OO 2 2 ^2 

<H.E[||T,(z.,)||^] + H _e-%— E[i.,((0,z)f]rfz+|M| -e-^^ ^E[(i.,((0, z])f ]dz 

/eo 2 2 ^2 /OO 2 2 ^2 

<H-E[||T,(^,)||g] + Hy _e-§._i^^z3+^dz + 4|Hy_ _e-§._dz 



£0 
00 2 



/OO 2 2 ^2 ; — 7 2 7^ 

<H.E[||T,(^,)||^] + Hy^ -^e-i^^KTz'+^dz + A\^,\ 1^ -e'^^ —dz 



eo 

2V2 

✓21 
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Finally we observe that for 77 > 



°° 2 1 1 2 

x^e"^ dx < -(77 + -)e"'' 
2 ?7 



and hence setting 77 — so that by assumption ?/ > 1 we have 



E[\-2fi{Ts{iys),d,T^{iys))o\] < ■E[\\Ts{u,)\\t] + C't^^ 



1x4 , '^2£0„_et 



J2 



where 



is a constant and 



23+1 



+ 00 



C2 = ImI 



8^2 



Lemma 3.8. -For 6 < Sq/2 we have 



1 — p V 0^ 



where C^,C2 and Eq are the same as in Lemma^ 
Proof. We have 



Jo 



00 / /"OO \ 

dxG's{x,y)vsidy)j dx 
d,,Gs{x,y) ==e"^^s^i^-— ^ ) i^sidy) ] dx 





00 / POO 



'2t:5 



dxGs{x,y)us{dy) ~ [ ^— e ^ ^ ^ i^s{dy) 



\ Jo V27rS d J 

2 (x+y)^ a; + 2/ 



:||9j;r5i/s||o - 2 / dxTsVsix) 



=e 2i 



Vs{dy)dx 



00 / /"OO 



,e * ^ Vs{dy) ) dx. 



/27r(5 



Moreover, 



^00 /'OO 
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p\ I Ts{vs){x)d{d^Ts{vs){x)) + j {d.,Ts{v,){x)Y dx 

•OO nOO 

= P[- I {d,Ts{v,){x)fdx + j {d,Ts{vs){x)f dx 
0. 



Also wc know 



(1 - p){Ts{v,),diTs{v,))^ = {l-p)[- {d.,Ts{vs){x)Ydx 



Therefore we have 



{Ts[vs),dlTs{u,))o+p\\d,n[vM 

2 (x+b)^ X + y 



-2p I dxTsVs{x) 

Jo Jo v27r(5 



e 2a 



Vs{dy)dx 



+P 



OO / poo 



< 



V^O 

OO 



V27r(5 / Jo 



2P / \/- — -dxTsUs{x) 
V P JO 



p 2 _ (x+«)2 X + y 

e 24 — — i>s{dy)dx 



OO / /'OO 



(x+y)2 X + y 



1 -p V2^ 

2 



\Jo 

OO 



/27r5 



= e 2a 



dx^il-p) {d^,Ts{vs){x)Ydx 







OO / />oo 



- ^/ — -dxT5Vs{x)] dx + p 

/o VV P / Jo \Jo y^-PV2TrS 



p 2 _ (x+y)2 a; + y 

e 2a 



OO / /"OO 



Vo 



/2^ 



^ iys{dy)j dx-{l-p) J {dxTs{vs)(x)) dx 



= (l-p) I {d^Tsv,{x)fdx+ ^ 



OO / /'OO 



1 - P io 

nOO 

(l-p) / (a,r5(i/,)(x))2dx 



==e ) 

v27rJ 



P 



OO / /"OO 



2 (x + y)2 X + ?/ / , X \ 

e 2a J c(a;. 

V 27rd 



l-P-Zo 

By the estimate for Pi obtained in Lemma 13.71 we have 

P 



1 — p 1 — pV 0^ 



P -L - P 

where , C2 and io are the same as in Lemma 13.71 

Now, combining Lemma 13.71 and 13.81 gives the following 
Theorem 3.9. Ifi>t is an M(R'^)-valued solution of iS.S}) and Zs(t) = T^vt, we have for 5 < 



nzsimi <\\zsm\i + h [nm^sWods + -^c'^sh 

Jo -L P 
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Corollary 3.10. Ifvt is a measure-valued solution of hS. 51) . then Vt G -ff'^, a.s. and IE||!/t||Q < oo, 
Vt > 0. 



Proof. By (|3.22p we have for small S that 



1 



nzsmi<\\zsm\\i + \i^\ / E\\Ts{,^s)\\ids , 

(1 - pW 

:=\\Zsm\l + \^i\ [ nZsis)\\lds + fiS,T), 



where 



/(<5,T) = - C'rS^T , 

1-p (l-p)')^ 



Applying Gronwall's inequality we have 

E\\Zsml < {\\Zsm\l + fiS,T))e\'^\'. 

It is clear that lini5_j.o f{S,T) ~ 0. Now let be a complete, orthonormal system for H'-' such 
that e Cfc(M+). Then by Fatou's lemma, 



E 





= E 







<liminfE||Z5(i)||2 < lliyolloe'^'', 

{5—^0 



Therefore i^t G and E||i/t||§ < oo, Vt > 0. 



□ 



Now we have proved the existence of an i^-density for the limit empirical measure vt , given 
that 1^0 has an L^-density. 



Theorem 3.11. Suppose that vq G iJ". Then \3. 5\) has at most one measure-valued solution. 

Proof. Let and be two measure- valued solutions with the same initial value vq , and both of 
them satisfy the boundary condition stated in Lemma [3.51 By Corollarv 13.101 vl^vf G a.s.. 
Let Vt = v\ — . Then Vt G and also Vt is a signed measure- valued solution to the equation 
p.5p . It is straightforward to extend all the estimates we have obtained to the case of the difference 
of two solutions as \vt\ < -\- i/f and the equations are linear. 

Therefore by the appropriate extension of Theorem 13.91 we have for 5 < e^/l 



2C2Teo 



-4/2S 
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As before, taking S ^ 0, we have 



E||i/t||2 < 1^1 f E\\\iys\\\lds ^ f E\\v,\\lds 
Jo Jo 



s I lo"^) 

and by Gronwall's inequality, we have = 0. □ 
This completes the proof of the uniqueness of the L^-valued solution to the equation p.Sp . 

3.4 The limit SPDE 

Substituting the Lebesgue representation for the empirical measure into (j3.5l) . integrating by 
parts and writing for the adjoint operator of A, we get 

(f>{x)v(t,x) dx = J (f>{x)v{0,x) dx + J J A(f>{x)v{s, x) dx ds 

+ J j ^(j)'{x)v{s,x) dxdMs 
t 



0(a:)w(O, x) + j J <j){x)A^v{s,x) dx ds 

~ lo I '^^^^^ ^^"'(■^'"^^^ dxdMs 
(l){x) (^v{0,x)+ J A^v{s,x)ds - J ■^{^v{s,x)) 



dM, dx 



As this holds Vcf) £ C we have shown that we have a weak solution to the SPDE given by 

/•* /•* g 

v(t,x) = v[U,x) + v{s,x)ds~ / — {y/pv{s,x)) dAIs, (3.23) 

Jo Jo dx 

with w(t, 0) = for all t e [0, T]. Alternatively, we can write this in differential form 

dv 1 d^v dv 

dv{t,x) = -fi—{t,x)dt + - — {t,x)dt - y^—{t,x)dMt, (3.24) 

with v{t, 0) = for all t £ [0, T] and v{0, x) = vo{x). This is a stochastic PDE that describes the 
evolution of the distance to default of an infinite portfolio of assets whose dynamics are given by 
(|1.2p . However the derivatives are only defined in the weak sense. 

We can now use the limiting empirical measure vt to approximate the loss distribution for 
a portfolio of fixed size N whose assets also follow (|1.2p . We do this by matching the initial 
conditions, thus setting 

1 ^ 



i=l 



where the Xq > 0, i = 1, . . . , N are the initial values for the distance to default of the assets in 
our fixed portfolio of size N. 
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3.5 Solving the SPDE 

The SPDE (|1.5p without the boundary condition is easily solved as 



v{t,x) =u{t,x- y/pMt), Va;eR,i>0, 



(3.26) 



where u(t, x) is the solution to the deterministic PDE 



'XX 




(3.27) 



with u(0,x) — vo{x). 

The SPDE with the boundary condition has been treated in [55]. This allows us to complete 
the proof of our existence and uniqueness theorem. 

Theorem 3.12. Let vo{x) G iJ^((0, oo)). The SPDE \1.5]) has a unique solution u e L'^{Q x 
{0,T),g,H\{0,oo))) and is such that xu^^ G L^{n x (0, T), ^^((o, oo))). 

Proof. The result follows from Theorem 2.1 of [28]. Thus all we have to do is ensure that the 
conditions of that Theorem hold in our setting. The boundary of the domain (0, oo) is the single 
point and hence we can take the function 'ip{x) min(a;, 1) in the Theorem. The single point 
boundary trivially satisfies the Hypothesis 2.1 of [35]. The coefficients of our SPDE are constants 
and hence satisfy the measurability requirement of Hypothesis 2.2 and the Lipschitz condition 
of Hypothesis 2.4. Hypothesis 2.3 also follows as the coefficients are constants and the initial 
condition is in i?^. □ 

Proof. ( of Theorem : Our previous work has shown that the empirical measure satisfies p.Sp 
and has a unique density in L^{{0, oo)). By Theorem 13 . 1 2 1 the SPDE with boundary condition has 
a unique solution in H^{{0, oo)). As this solution satisfies (|3.5p . by the uniqueness of solutions, it 
must be the density for our empirical measure. Thus our density satisfies the SPDE. □ 

We note that we can derive a formal expression for Lt in terms of the density after integrating 
by parts. 



Lt =1 



1 



L 
L 



v(t, x)dx 



( 



t;(0,x)- / 



Jo 




Vxx{s,x)ds 



1 




ds 
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Since > 0,Vi and XI is a continuous process, we can conclude that Tq > 0,Vi. Thus 

1 ^ 

i=l 

therefore 

f +CXD 



D{R+ U {0}) 1 = / w(0, x)da;. 
Jo 

Moreover we have v{s,x) — >■ 0,Vx{s,x) — > 0, as a; — > cxo and w(s, 0) = 0, Vs. Therefore, provided 
that Vx{s, 0), the right derivative of v{s, x) with respect to x at the point a; = 0, exists we would 
have ^ 

Lt J Vx{s,0)ds. 

One issue that has not been addressed is the existence of solutions to this equation. We 
note that the work of Lototsky [3^ shows that there is a classical solution to this SPDE over a 
bounded domain (0, K), with Dirichlet boundary conditions at and K, provided that the initial 
condition is smooth enough. 



3.6 The portfolio loss 

We would like to price portfolio credit derivatives whose values depend on the cumulative 
defaults occurring within a reference basket of risky assets. The key to pricing these instruments 
is determining the joint loss distribution. We have just derived an equation that describes the 
evolution of the empirical measure of the limiting large portfolio of assets. At any future value 
in time, we can determine the loss in the portfolio by calculating the total mass of the empirical 
measure of assets that have not defaulted. Thus the portfolio loss can be approximated by 

Lf = NLt, 

where N is the number of assets in the portfolio. We note that given the initial condition p.25p 
we have Lq = 0. Also, due to the way in which defaults are incorporated into the model, we have 

< it < 1, for t > 

P{Ls >K)< P{Lt >K), for s < t, 

which ensures that there is no arbitrage in the loss distribution. Both of these properties are 
expected for a model of cumulative loss in a portfolio. 



3.7 A connection with filtering 

We note that the SPDE can be viewed as a PDE with a Brownian drift. This is easily 
seen through an interpretation as the Zakai equation for a filtering problem. Let {il, J-, P) be a 
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probability space. Under P we define the signal process X to be a stochastic process satisfying 



dX = fidt - y^dM + ^1 - pdW, t < tq 
Xt = 0, t>To 

where tq = mi{t : Xt = 0}, where fx, p are constants and M and W are independent Brownian 
motions and Xq — x. The observation process Y is taken to be just the market noise, 

dYt - dMt, 

then the Zakai equation (see for example [2]) for the conditional distribution of the signal given 
the observations is exactly our SPDE. 

Thus, by standard filtering theory, if we want to compute a functional of the signal we need 
to calculate 

m^(t) =E(^/.(Xt)|-Ff ) = j ^{y)u{t,x)dx. 

This means that the probability distribution for the position of a company given the market noise 
has a density u(t,x) satisfying 

du{t,x) = {-p,Ux{t,x) + ^Uxx{t,x))dt - ^Ux{t,x)dMt, 

with m(0, x) = uo{x), that is the initial guess at Xq is the density uo{x) and u{t, 0) — 0. Thus for 
the loss function we are interested in computing the proportion of companies that have defaulted 
by time t and this can be found by computing m^(t) for iplt) — /{ro<t}- If we start from a given 
fixed point so that ^0(2:) is a delta function at x. Then 

Lt = mj,{t) = F(inf Xs < 0| J"f ). 

s<t 

Now the process X can be written as a Brownian motion with drift 

Xt^x + fit- y^Mt + ^1-pWt, 
and if we are given M , this can be expressed as 

where f{t) = fit — y/pMt is a deterministic time dependent drift function which is a fixed random 
path. 

Thus to compute the random loss function we set x' = — p, g{t) = f{t)/^J\ — p and 
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write 



P"(inf X, < 0|J-f ) 

s<t 



P (inf x' + g{s) + Ws < 0| J"f ) 



's<t 



Piinf g{s) + Ws < -x'\Tf^). 

S<t 



In the case where we have a general initial distribution uo{x)^ the loss function is then 



/ 

^0 



■oo 



Lt = 



Mo(x)P(inf gs + Ws < -x/y/1 - p\Tf^)dx. 



's<t 



Thus wc can try to compute this by solving the hitting time problem for Brownian motion 
with time dependent drift for a fixed realization of the market noise. It is straightforward to use 
this to simulate a realization of the loss function. 

To derive this SPDE we made some simplifying assumptions. The first of these arose when 
specifying the asset processes in ()1.2|) . We had to set the drift and volatility of all the assets to 
some common value. For the drift this is not a problem, because under the risk neutral measure it 
will be transformed to a value that excludes arbitrage. The fact that there is only one yield curve 
means that this value will be the same for all assets. If our reference portfolio contained entities 
denominated in more than one currency this would not be the case and some approximation would 
have to be made. 

This argument cannot be used for the volatility as it is not affected by a change of measure. 
Therefore, it would seem that giving the assets one common value of volatility is a very restrictive 
assumption. However, for any given value of the volatility we still have the freedom to choose the 
default barrier specific to any one asset. Via the distance-to-default transformation this freedom 
manifests itself in our particular choice of starting value for each process. The effect of changing 
the barrier and changing the volatility is very similar. To see this note that default risk is measured 
by how many standard deviations away from the barrier our process is. To increase the default 
risk we need to reduce this distance which can be done by either increasing the standard deviation 
or moving the barrier closer. Although these are clearly not equivalent transformations they have 
a very similar effect and so the single volatility assumption is not as restrictive as it initially 
appears. 

Having a single volatility number also eases calibration as we do not have to estimate the 
volatilities of all of the entities within our portfolio. Instead, we will have to replace it by some 
'average' market volatility. Not only will this help day-to-day calibration stability but it means 
that credit derivative prices will be a function of one volatility parameter only. This is usually 
a desirable property from a practitioner's point of view as it allows one to take a view on that 
parameter; this cannot be done if there were a single parameter for each entity within our portfolio. 

The major simplification that allowed us to derive our SPDE came when we moved to an 
infinite dimensional limit. In this limit, the idiosyncratic noise of the individual assets is averaged 
out. In fact, we could have any number of idiosyncratic components, provided they are indepen- 
dent and uncorrelated, and they would average out to zero. It is only the correlated components 
between the assets that remain i.e. the market risk. Note that this means that if the limiting 
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portfolio was fully diversified, that is had no correlation, there would be no noise in the limit and 
the limit portfolio would evolve deterministically! 



4 Numerical solution 

We outline in the following a numerical method for approximating the solution to the SPDE, 
which we use in the market pricing examples in the next section. We start with the SPDE (|3.5|) 
in weak form, repeated here for convenience, 

((/), Vt) = (0, vq) + \ (^0, Vs) ds+ f (y/pcl)', Vs) dMs 
Jo Jo 

for almost all t and all smooth test functions <j> ^ C. It follows from Theorem II . II that vt has as 
one component the density v (describing the non-absorbed element) satisfying 

(0,7;(t,-)) = (0,«(O,-))+ / {A(p,v{s,-))ds + ^p [ (0',z;(s,.))dM„ (4.1) 

Jo Jo 

where here we write (•,•) for the inner product. Integrating by parts, noting from Theorem 
13. 121 that v{t, ■) e Hi with dense subspace (7, 

((/),«(*,•))+/ a((/),«(s,-))ds = (0,i'(O,-)) + VP / {<i}'Ms,-))dM, 
Jo Jo 



for all (j) £ i?o, where 



a{<l,,v) = ]^{<j,',v')-^p{cj,',v). 



4.1 Finite element approximation 

Let Vh C HI{[xq,xn]) be the space of piecewise linear functions on a grid Xi < . . . < xn, 
which are zero at xi = and a; at a sufficiently large value fsee 14. 3|) . Denote further by {0„ : 
1 < rt < iV} the standard finite element basis (see e.g. [37] for standard finite element theory and 
approximations to PDEs). Restricting both the solution and test functions to V^, 

{4>n,Vh{t,-)) + a{(t)n,Vh{s,-))ds = {(t)n,Vh{0,-)) + ^ {(j)'^, Vh{s , ■)) dM s 

Jo Jo 

(for all 1 < n < N) defines a semi-discrete finite element approximation. 

Using the stochastic ^-scheme (see [18]) for the time discretisation of the resulting SDE system, 

{<l^n,v^+^) + 6 At a(0„, 1;^+') = O - (1 - 0)At a(0„, <) + (0'„, <) TAtl*™, (4.2) 

where $m~^(0, 1), Ai = tm+i — tm is assumed constant and — X^^lLi ^nVn- Thus one gets 
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a linear system 

{M + eAtA)v"'+^ = [M - (1 - e)MA)v'^ + 7p\/Ai$„L>w", (4.3) 
where = (u™, . . . , w]^) and the standard finite element matrices are given by 

A^j = a{(j>i,(j)j), l<i,j<N, 
A, = (<^-,</'j), l<»,j<^. 

This gives a pathwise (in M , the market factor) approximation to the SPDE solution via timestep- 
ping from an initial density Vh{0, •), which is found by projection of VNf.t from ()1.4p with Nf 
firms onto the finite element space (see e.g. [36], [38]). 



4.2 Simulating tranche spreads 

For a given (numerical) realisation of the market factor, we can approximate the loss functional 
LTfc at time Tk by 

fXM N-i 

ir, = 1 " / Vh{Tk,x) dx^l~hY,vZ' (4.4) 
■^0 n=l 

where m = T^/At. If we explicitly include the dependence on the Monte Carlo samples $ — 

{^i)i<i<i in -^Tfc(*^)' where $i as in (|4.2p are drawn independently from a standard normal 

distribution, then for Ngims simulations with samples = {^i)i<i<i , 1 < ^ < Nsimsi we simulate 

the outstanding tranche notional (|2.3I) as 

[Zt, ] w £"0 [max(d - L^^ , 0) - max(a - Li^^ , 0)] 

« 5] (max(d-L$,J$'),0)-max(a-L^^(<I>'),0)). 

sims ^ ^ 

This gives simulated tranche spreads via (|2.5p . (|2.6I) and (|2.7I) . 



4.3 Accuracy and further approximations 

We now discuss the approximations made previously and further simplifications made in the 
numerical implementation of the examples in the next section. 

It is necessary for the finite element discretisation to approximate the semi-infinite boundary 
value problem for the SPDE by one on a finite domain. It is expected that if the upper bound- 
ary is sufficiently large, dependent on the initial distances-to-default and model parameters, the 
probability of crossing this boundary can be made negligible and zero boundary conditions are 
appropriate. We have checked this to be the case for the following numerical simulations but do 
not have a theoretical justification at this point. 

The derivation of the SPDE and finite element solution assume initial data, however in 
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practice we want to use a sum of atomic measures (|3.25p corresponding to the distance-to-default 
of individual firms, as backed out from CDS spreads. We deal with this by projecting these data 
onto the finite element basis (see e.g. [3^, [55]). 

The majority of the literature on stochastic finite element methods deals with stochastic 
diffusion coefficients (see e.g. [9] and subsequent work) and we are not aware of results which 
cover our setting with stochastic drift. From standard finite element approximation results for 
PDEs (see e.g. [37]), one would expect (pathwise) convergence order two in h for solutions in H^, 
but Theorem 13.121 suggests weaker regularity at the absorbing boundary, which we also observe 
in the numerical solutions. This does not show a measurable impact on the numerical accuracy 
in practice. The weak approximation order of the Euler scheme for SDEs, and that for the chosen 
fully implicit scheme for PDEs (6* = 1 in (|4.2p ). is one (in At). In this case, the scheme is stable 
in the mean-square sense of [18| . This is confirmed by numerical experiments, but a rigorous 
numerical analysis is beyond the scope of this paper. 

A common approximation to the finite element system is to 'lump' M in (|4.3p in diagonal form, 
interpretable as application of a quadrature rule, and ultimately results in M being replaced by 
a multiple of the identity matrix. With this approximation, the finite element scheme becomes 
identical to a central finite difference approximation. 

A further simplification is suggested by the solution p.26p of the SPDE without absorbing 
boundary condition, which decouples the solution into the PDE solution (|3.27p on a doubly- 
infinite domain, and a random (normal) offset. This is easy to implement if we apply boundary 
conditions only at a discrete set of times. In analogy to discretely sampled barrier options, this 
corresponds to a situation where we observe default not continuously, but only at discrete dates. 
The numerical results in the next section were obtained in this way with default monitoring at 
payment dates for computational convenience. This introduces a small shift in the calibrated 
parameters compared to the SPDE with continuously absorbing barrier but the reported results 
on tranche spreads are almost identical. 

The Monte Carlo estimates of outstanding tranche notionals and subsequently tranche spreads 

— 1/2 

converge per N^^J^^ . The variance relative to the spread is larger for senior tranches due to the 
rarity of losses in these tranches, as illustrated by Figure |T] Importance sampling could cure this 
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Figure 1: Monte Carlo estimators with standard error bars for expected losses (|2.4|) in tranches 
[0, 3%], [6%, 9%], [12%, 22%], for Ns^ms = 16 • 4'^-!, fc = 1, 10, and a typical set of parameters, 
maturity T = 5. 

problem but was not found necessary for the purposes of this study. 
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Numerical parameters were in the following adjusted such that the (heuristically) estimated 
approximation error was sufRciently small compared to the effects observed by varying model 
parameters. 

5 Market pricing examples 

5.1 Calibration to index tranches 

In this section, we analyse our model's ability to price regular index tranches for all maturities 
and investigate the implied correlation skew. We consider performance pre and post the onset 
of the credit crunch, illustrating the model's inherent ability to cope with a variety of credit 
environments. 

Throughout the analysis, we infer the initial condition from market spreads for the underlying 
index constituents, rather than allowing it to be a free parameter to be fixed by calibration to 
index tranches. This is to be consistent with CDS spreads for the individual constituents. We 
do this by backing out the distance-to-default for each constituent from its five-year CDS spread 
and then aggregating these. Note that as we model the distance-to-default as in (|l.ip . different 
volatilities of the underlying firms can be taken into account by rescaling. As a consequence, the 
initial condition is driven by both the level of constituent spreads and their dispersion. 

We study the ability of our model to price index tranches on two dates: February 22, 2007 and 
December 5, 2008. These dates are chosen specifically to investigate the flexibility of the model 
to cope with different market and spread environments. February 22, 2007 was pre-crisis when 
spreads were tight and curves upward sloping; December 5, 2008 was at the height of market 
volatility, when spreads were at their widest and curves frequently inverted. 

We set R = 40%, the level typically assumed by the market for investment grade names, and 
for each date, calibrate the model to 5, 7 and 10-year index spreads using the volatility, a. r is 
the risk-free rate obtained from the Euro swap curve. (N.B. the correlation parameter, p, does 
not come into this calibration since index spreads depend only on the expected losses, which are 
identical to the sum of default probabilities and hence correlation- independent.) 

Table [1] shows the traded and model index spreads for Feb 22, 2007. Since we derive the 
initial condition from constituent spreads, we only have one free parameter, the volatility a, for 
calibrating all three index spreads. Increasing a to increase model spreads also causes the initial 
distance-to-default for each constituent to increase (since CDS spreads are fixed), so index and 
tranche spreads are less sensitive to changes in volatility than they would be if the initial condition 
was specified independently. 

Table [3] shows the same results for Dec 5, 2008. In this highly distressed state, we notice that 
spreads are dramatically wider and the curve is inverted with 5-year > 7-year > 10-year spreads. 
Our simple model again does a good job of calibrating all three spreads. This is achieved by a 
smaller distance-to-default for the initial positions in combination with a lower volatility, triggering 
more defaults in the near future. The 5-year point is a little low, which is a shortcoming of using 
a purely diffusive driving process: it can be hard to generate sufficient short-term losses. We refer 
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Maturity Date 


Fixed Coupon (bp) 


Traded Spread (bp) 


Model Spread (bp) 


20/12/2011 


30 


21 


19.6 


20/12/2013 


40 


30 


30.7 


20/12/2016 


50 


41 


41.0 



Table 1: The fixed coupons, traded spreads and model spreads for the iTraxx Main Series 6 index 
on February 22, 2007. Parameters used for the model spreads are r = 0.042, a = 0.22, R = 0.4. 

to Section |6] for a discussion of extensions to jump and stochastic volatility driven processes. 

For the parameters from the calibration in Table[l] Tabled illustrates the correlation sensitivity 
of the 5, 7 and 10-year index tranches in the pre-crunch environment. We note that model spreads 
illustrate the behaviour we would anticipate: 

• Equity tranche spreads decline with increasing correlation whilst spreads for other tranches 
generally increase with correlation. As correlation increases, there are less likely to be a 
few defaults, and so the equity tranche becomes less risky and its spread decreases. The 
probability of a greater number of defaults increases with increasing correlation and so 
spreads on the more senior tranches increase with correlation. 

• A notable exception is the 10-year junior mezzanine tranche (3% — 6%) which behaves more 
like an equity tranche and has declining spreads with increasing correlation. This is because, 
for the parameters used, the expected index loss is between 3% and 6%. The risk of this 
tranche therefore decreases, along with the spread, as correlation increases, making losses 
in this tranche less likely. 

• The 7-year junior mezzanine tranche (3% — 6%) spreads indicate the transition, as maturity 
increases, from positive to negative correlation sensitivity by exhibiting a humped shape. 

• For the 5 and 7- year junior mezzanine and 10-year senior mezzanine tranches, spreads decline 
with increasing correlation for high values of correlation. 

Figure [2] illustrates the 5, 7 and 10-year implied correlation skew - the value of correlation 
that gives a model spread equal to the market spread for each tranche and maturity. 

• With the exception of the 0% — 3% tranche, we see similar behaviour and levels for all 
three maturities. This consistency across the term-structure suggests that the dynamics 
underlying the model are realistic, even in its simple form. 

• 5-year implied correlations are generally high relative to the others and 10-year values rela- 
tively low. To achieve consistency of the correlation parameter across maturities, a driving 
process with the ability to generate more default events in the short-term would be required, 
eg a more general Levy process for the market factor. 

• An anomaly is revealed by the 3% — 6% implied correlations and the corresponding row data 
in Table [21 where it is seen that the correlation dependence of model tranche spreads flips 
from increasing to hump-shaped to decreasing for maturities running from 5 to 10 years. This 
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5 Year 


Tranche 


Market 


p = 0.1 


p = 0.2 


p = 0.3 


p = 0.4 


p = 0.5 


p = 0.6 


p = 0.7 


0%-3% 


7.19 % 


7.55 % 


4.99 % 


2.14 % 


-0.71 % 


-3.48 % 


-6.17 % 


-8.78 % 


3%-6% 


41 


15.6 


55.6 


86.4 


106.1 


116.2 


119.5 


117.4 


6%-9% 


10.8 


0.7 


9.1 


25 


40.3 


54.5 


65.2 


71.7 


9%-12% 


5 





2.2 


8.2 


18.8 


28.6 


37.2 


45.4 


12%-22% 


1.8 





0.2 


1.7 


4.9 


9.8 


16.1 


22.5 


22%-100% 


0.9 











0.1 


0.3 


0.7 


1.5 




7 Year 


Tranche 


Market 


p = 0.1 


p = 0.2 


p = 0.3 


p = 0.4 


p = 0.5 


p = 0.6 


p = 0.7 


0%-3% 


22.1 % 


27.45 % 


19.97 % 


13.79 % 


8.31 % 


3.27 % 


-1.47 % 


-6.04 % 


3%-6% 


110 


130.6 


183.3 


202.2 


206 


201.5 


191.6 


177.8 


6%-9% 


32.5 


15.3 


52.4 


80.5 


99.1 


110.6 


116.1 


116.9 


9%-12% 


15 


1.8 


17.4 


37.1 


54.3 


67.1 


76.5 


82.7 


12%-22% 


4.9 


0.1 


2.3 


8.9 


19 


29.9 


39.5 


47.9 


22%-100% 


2 








0.1 


0.4 


1.1 


2.3 


4.1 




10 Year 


Tranche 


Market 


p = 0.1 


p = 0.2 


p = 0.3 


p = 0.4 


p = 0.5 


p = 0.6 


p = 0.7 


0%-3% 


38 % 


42.51 % 


32.51 % 


24.13 % 


16.65 % 


9.71 % 


3.11 % 


-3.33 % 


3%-6% 


302.5 


375.8 


354.9 


331.9 


308.1 


283.5 


258.3 


231.9 


6%-9% 


83 


101.4 


147.3 


166.2 


173.6 


174.4 


170.8 


163.8 


9%-12% 


37 


24.3 


64.1 


90.9 


107.7 


117.8 


122.9 


124.1 


12%-22% 


12.5 


2 


13.5 


29.1 


44.4 


57.5 


68.2 


76.5 


22%-100% 


3.6 





0.1 


0.6 


1.5 


3 


5.1 


7.7 



Table 2: Model tranche spreads (bp) for varying values of the correlation parameter. The equity 
tranches are quoted as an upfront assuming a 500bp running spread. The model is calibrated to 
the iTraxx Main Series 6 index for Feb 22, 2007. Market levels shown are for this date; model 
parameters are r = 0.042, a = 0.22, R = 0.4. 

has the following effect: for 5 years, there is a unique implied correlation for this tranche; 
for 7 years, a second, higher, correlation (just under 1) also fits this tranche; for 10 years, 
only a single high correlation can fit the market spread. Essentially, the implied correlation 
curves in Figure [2] are shifted downwards with increasing maturity. The alternative higher 
branches, where applicable, are not included in the Figure. When a curve crosses zero (in 
the case of the 10-year 3% — 6% tranche), we have set the implied correlation to zero (instead 
of the value of around 0.42 from the higher branch which exactly reproduces the market 
quote). For pricing and (especially) hedging purposes, continuous dependence of implied 
correlations with respect to maturity and market data is clearly desirable. The lack of a 
calibration which is both stable and exact underlines the need for a richer model. 

Table |4] shows the correlation sensitivity of the Dec 5, 2008 index tranches with parameters 
from the calibration in Table [S] We notice that relative to Table [21 spreads are highly distressed, 
the index is inverted and tranche spreads are flat to inverted across maturities. As a result, 
the tranches exhibit very different sensitivity to correlation than before, however there are some 
common themes and extensions to earlier behaviour; 
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Figure 2: Implied Correlation Skew for iTraxx Main Series 6 Tranches, Feb 22, 2007. 




The implied correlation for each tranche is the value of correlation that gives a model tranche spread 
equal to the market tranche spread given in Table [2] Model parameters are r = 0.042, a = 0.22, R — 0.4. 



Maturity Date 


Fixed Coupon (bp) 


Traded Spread (bp) 


Model Spread (bp) 


20/12/2013 


120 


215 


207 


20/12/2015 


125 


195 


195 


20/12/2018 


130 


175 


176 



Table 3: The fixed coupons, traded spreads and model spreads for the iTraxx Main Series 10 
index on December 5, 2008. Parameters used for the model spreads are r = 0.033, a ~ 0.136, 
R = 0.4. 



• Default probabilities for the index and its constituents are very high. The index expected 
loss is therefore much greater than before, illustrated by the fact the first three 5-year 
tranches and the first four 7 and 10-year tranches have declining spreads with increasing 
correlation. This contrasts with just the equity and 10-year junior mezzanine tranches in 
Feb 2007. 

• Much higher levels of p are needed to replicate market prices than in pre-crunch times, 
consistent with the fact that systematic risk is a much greater concern at this time. 

• Too much of our model's portfolio loss distribution lies in the middle tranches: 6% — 22%; 
more weight needs to be in the tail to be able to replicate 22% — 100% tranche values. 
The same model shortcoming holds for all maturities and reflects the need for a more 
sophisticated driving process. 
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5 Year 


Tranche 


Market 


p = 0.3 


p = 0.4 


p = 0.5 


p = 0.6 


p = 0.7 


p = 0.8 


p = 0.9 


0%-3% 


71.5 % 


81.88 % 


75.9 % 


69.56 % 


63.02 % 


56.25 % 


49.16 % 


41.65 % 


3%-6% 


1576.3 


2275.2 


1978.5 


1743.2 


1546.8 


1374.6 


1222.8 


1090.1 


6%-9% 


811.5 


1273.1 


1168.2 


1079.7 


1001.4 


931.3 


864.6 


796.3 


9%-12% 


506.1 


775.7 


765.8 


748.6 


724.7 


695.8 


663.2 


629.1 


12%-22% 


180.3 


307.8 


353.3 


384.7 


405.5 


418.1 


423.4 


420.5 


22%- 100% 


77.9 


9.2 


16.5 


25 


34.3 


44.5 


55.7 


68.1 




7 Year 


Tranche 


Market 


p = 0.3 


p = 0.4 


p = 0.5 


p = 0.6 


p = 0.7 


p = 0.8 


p= 0.9 


0%-3% 


72.9 % 


84.03 % 


78.98 % 


73.26 % 


66.93 % 


60 % 


52.41 % 


44.13 % 


3%-6% 


1473.2 


2327.3 


1985.7 


1715.2 


1493.4 


1308 


1147.8 


1001.3 


6%-9% 


804.2 


1344.2 


1199 


1085.2 


988.2 


900.7 


820.9 


747.9 


9%-12% 


512.4 


855.4 


808.4 


765.3 


725.3 


684.8 


643 


600.4 


12%-22% 


182.6 


375.4 


401.7 


417.6 


425.6 


427.4 


423.1 


411.8 


22%- 100% 


75.8 


14 


22 


30.6 


39.6 


49.3 


59.7 


71.2 




10 Year 


Tranche 


Market 


p = 0.3 


p = 0.4 


p = 0.5 


p = 0.6 


p = 0.7 


p = 0.8 


p = 0.9 


0%-3% 


73.8 % 


85.13 % 


80.57 % 


74.99 % 


68.51 % 


61.31 % 


53.31 % 


44.22 % 


3%-6% 


1385.5 


2270.8 


1895.7 


1611.1 


1385.8 


1195.3 


1032 


889.6 


6%-9% 


824.7 


1332.2 


1164.2 


1033.7 


925.5 


833.5 


749.8 


669.7 


9%-12% 


526.1 


870.8 


798.8 


740.7 


689.3 


640.5 


592.1 


543.1 


12%-22% 


174.1 


406.1 


414.9 


417.5 


415.6 


409.8 


400.2 


385.3 


22%- 100% 


76.3 


18.3 


26.1 


34 


42.1 


50.6 


59.7 


69.8 



Table 4: Model tranche spreads (bp) for varying values of the correlation parameter. The equity 
tranches are quoted as an upfront assuming a 500bp running spread. The model is calibrated to 
the iTraxx Main Scries 10 index for Dec 5, 2008. Market levels shown are for this date; model 
parameters are r = 0.033, a = 0.136, R = 0.4. 



5.2 Forward steirting CDO contracts 

These contracts are obligations to buy or sell protection on a specified tranche for a specified 
spread at some specified time in the future. Although these instruments are traded infrequently, 
their pricing and hedging is an active research topic. We will look at two types of forward starting 
CDO: one that resets the cumulative loss at the forward start date and one that docs not. For 
discussion purposes we will refer to these as resetting and non-resetting respectively but it should 
be borne in mind that these are not standard market terms. 

5.2.1 Non- resetting forward CDO tranche 

For a non- resetting forward CDO tranche defined over the time interval [T, T*], the cumulative 
losses incurred up to time T count towards the total loss in the tranche for all t with T <t. This 
feature makes pricing straightforward and analogous to a forward CDS contract. 

Consider a portfolio with m entities in the reference portfolio. We define the total loss on the 
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portfolio at time t by 

m 

T.^^Mn<t}- (5.1) 



i=l 



If the forward tranche has attachment point a and detachment point d then the outstanding 
tranche notional, Zt, is given as 

Zt = [d - Lt]+ - [a - Lt] + . (5.2) 

The value of the forward tranche contract is again given by the difference between the fee leg and 
the protection leg. So far the setup has been the same as the standard CDO tranche. The only 
difference when pricing this forward contract is the fact that now we are only interested in the 
payment dates T^, i = 1, . . . , n where T < Ti <...,< r„ < T*. Using these payment dates the 
present value of the coupon payments given a forward spread s is 

«^'^^ = ^E6|^IEQ[ZtJ. (5.3) 

The protection leg is given by 

Today, the value of the forward starting contract is zero and hence the forward break-even spread 
is given by 

s = —^. (5.5) 

yfee ^ ' 

5.2.2 Resetting forward CDO tranche 

With this contract, the cumulative loss up to time T is ignored and the value of the tranche is 
dependent only on the further loss incurred after time T. If the forward tranche has attachment 
a and detachment d then this is equivalent to a non-resetting forward tranche with attachment 
{Lt + a) and detachment (Lt + d) . Using the same payment dates as the non-resetting forward 
contract we define the effective forward loss at time Ti by 

= — Lt, (5.6) 

which gives the forward tranche notional as 

Zt^[d-Lt]+ -[a-Lt] + . (5.7) 



With these new definitions the forward break-even spread can be calculated as before using (|5.3I) 
(Ea and (1531). 
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5.3 Forward pricing results for the pre-crunch state early 2007 

We value resetting and non-resetting forward CDO contracts for a range of forward starting 
dates T. The data used is for the European iTraxx Main Series 6 index from February 22 2007. 
The index fixed coupons and traded spreads are shown in table [1] and we use a constant risk- free 
rate of 4.2% obtained from the Euro swap curve. The tenor of the forward contracts is always 
five years i.e. T* — T = 5. The forward dates we use are years i.e. the spot spread and the 1, 
3 and 5 year forward starting dates. The forward break-even spreads for the non-resetting and 
resetting forwards are shown in table [5l 



p = 0.1 


non-resetting 


resetting 




T = 


T = 1 


r = 3 


T = 5 


T = 1 


T = 3 


r = 5 


0%-3% 


7.55 % 


24.06 % 


51.32 % 


64.31 % 


24.05 % 


50.26 % 


60.95 % 


3%-6% 


15.6 


71.0 


374.2 


911 


70.9 


312.3 


484.0 


6%-9% 


0.7 


4.7 


58.7 


231 


4.7 


43.5 


79.5 


9%-12% 





0.7 


9.7 


54 


0.7 


6.6 


10.3 


12%-22% 








0.7 


4.5 





0.3 


0.5 


22%-100% 























p = 0.3 






0%-3% 


2.14 % 


13.87 % 


31.87 % 


36.77 % 


13.87 % 


33.04 % 


42.44 % 


3%-6% 


86.4 


174.2 


441.9 


704.3 


174.1 


402.9 


545.3 


6%-9% 


25.0 


62.6 


188.5 


358.6 


62.6 


164.1 


231.1 


9%-12% 


8.2 


24.1 


88.1 


197.6 


24.1 


74.4 


107.9 


12%-22% 


1.7 


5.0 


24.3 


63.2 


5.0 


19.9 


25.8 


22%-100% 





0.1 


0.4 


1.3 





0.2 


0.2 


p = 0.5 






0%-3% 


-3.48 % 


5.18 % 


17.96 % 


20.16 % 


5.19 % 


20.11 % 


28.31 % 


3%-6% 


116.2 


193.9 


400.8 


532.2 


193.8 


384.5 


507.8 


6%-9% 


54.5 


101.6 


228.2 


341.4 


101.5 


211.2 


283.7 


9%-12% 


28.6 


59.8 


142.9 


237.0 


59.7 


128.9 


177.0 


12%-22% 


9.8 


22.0 


62.0 


118.5 


22.0 


54.4 


73.9 


22%-100% 


0.3 


0.8 


2.9 


6.4 


0.8 


2.3 


2.1 



Table 5: The non-resetting and resetting forward spreads (bp) for varying values of the correlation 
parameter. The equity tranches are quoted as an upfront assuming a 500bp running spread. The 
model is calibrated to the iTraxx Main Series 6 index for 22 Feb 2007. All forwards have a tenor 
of 5 years. 



5.3.1 Non-resetting forward CDO tranche 

First we focus on the non- resetting tranches. From table [5] we observe the following points: 

• As the forward start date increases, the break-even forward spread increases for all tranches. 

• For the junior mezzanine (3-6%) tranche, as the forward start date increases the spread 
sensitivity to correlation changes sign. The sensitivities of all other tranches are single 
signed. 
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Figure 3: Rate of expected losses in tranehes [0,3%], [3%,6%], [6%,9%], [9%, 12%] in 2007, for 
correlation p — 0.5, for non-resetting (left) and resetting (right) losses. 

Both of these observations can be explained by the fact that losses in the portfolio are cumulative. 
As time passes, the total loss in the portfolio accumulates and so the attachment and detachment 
points of non-resetting forward tranches effectively move down the capital structure. In other 
words, forward equity tranches start behaving like very narrow spot equity tranches, forward 
junior mezzanine tranches start behaving like spot equity tranches and so on. On the forward 
start date, investors will require additional compensation for holding these now riskier tranches 
and so the break-even forward spread increases. 

This also gives the reason why the correlation sensitivity of the junior mezzanine tranche 
changes sign. For a start date sufficiently far into the future, the tranche is expected to be an 
equity tranche which has a negative correlation sensitivity. 

5.3.2 Resetting forvirard CDO tranche 

Now we turn to the resetting tranches which from a dynamic modelling point of view can be 
considered the more interesting of the forward contracts. From table [5] we observe the following 
points: 

• For the forward tranches with a start date in 1 years time the break-even spreads are the 
same as for the non-resetting forward tranches. 

• As the forward start date increases, the break-even forward spread date generally increases. 

Addressing these observations in order, the reason for the first point is the nature of the structural 
model. Because of the diffusive nature of the asset processes, the probability of defaults occurring 
in the short term is very low. In Section 3 this was highlighted as one of the major downsides for 
this type of model. The consequence of that property here is that the cumulative loss within the 
first year is negligible and so we have Lt « Lt . This in turn leads to the same break-even spreads 
for both types of forward contract. 

The second point can be explained simply by the potential for a decrease in credit quality due 
to the natural diffusion of the asset processes. This is also present in the non-resetting tranche 
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prices but there the increase in spreads is dominated by the move down the capital structure. An 
exception is the super senior (22-100%) tranche for high correlation. We come back to this when 
discussing the distressed state where this effect is more pronounced. 

5.4 Forward pricing results for the distressed state late 2008 

We again value resetting and non-resetting forward CDO contracts for a range of forward 
starting dates T . The data used is for the European iTraxx Main Series 6 index from December 
5 2008. The index fixed coupons and traded spreads are shown in table [3] and we use a constant 
risk-free rate of 3.3% obtained from the Euro swap curve. The contract details are as before. 
Forward break-even spreads for the non-resetting and resetting forwards are shown in tables [H 



p = 0.3 


non-resetting 


resetting 




r = 


r = 1 


T = 3 


T = 5 


r = 1 


r = 3 


T = 5 


0%-3% 


81.88 % 


82.40 % 


69.26 % 


56.96 % 


83.21 % 


77.69 % 


68.98 % 


3%-6% 


2275.2 


3106.5 


2658.7 


1989.3 


2479.9 


1830.2 


1206.6 


6%-9% 


1273.1 


1708.0 


1771.4 


1401.2 


1348.2 


901.9 


522.1 


9%-12% 


775.7 


1045.9 


1205.5 


1024.9 


811.7 


487.7 


238.9 


12%-22% 


307.8 


425.3 


570.1 


546.6 


313.3 


143.9 


51.4 


22%-100% 


9.2 


13.9 


24.1 


29.4 


8.0 


1.8 


0.2 


p = 0.5 






0%-3% 


69.56 % 


66.22 % 


47.22 % 


34.18 % 


70.09 % 


63.07 % 


53.86 % 


3%-6% 


1743.2 


2090.3 


1580.0 


1119.3 


1858.6 


1441.5 


1033.5 


6%-9% 


1079.7 


1341.4 


1149.0 


869.2 


1141.9 


845.4 


565.4 


9%-12% 


748.6 


938.6 


890.4 


690.9 


787.1 


549.8 


333.5 


12%-22% 


384.7 


495.5 


532.2 


455.6 


400.0 


241.7 


123.6 


22%-100% 


25 


33.4 


44.8 


44.8 


23.3 


8.6 


2.5 


p = 0.7 






0%-3% 


56.25 % 


49.41 % 


29.70 % 


16.99 % 


55.67 % 


47.50 % 


38.76 % 


3%-6% 


1374.6 


1481.6 


1040.3 


734.7 


1428.3 


1146.6 


864.3 


6%-9% 


931.3 


1072.1 


827.8 


602.8 


973.7 


766.2 


549.2 


9%-12% 


695.8 


818.8 


689.6 


516.9 


731.3 


554.3 


379.7 


12%-22% 


418.1 


513.8 


478.1 


384.8 


443.0 


310.4 


188.1 


22%-100% 


44.5 


55.9 


64.1 


57.9 


42.7 


20.8 


9.5 



Table 6: The non-resetting and resetting forward spreads (bp) for varying values of the correlation 
parameter. The equity tranches are quoted as an upfront assuming a 500bp running spread. The 
model is calibrated to the iTraxx Main Series 6 index for 5 Dec 2008. All forwards have a tenor 
of 5 years. 

Comparing figures S] with [31 the expected loss rate peaks at the short end in the distressed 
environment of 2008. The effect on the resetting forward starting CDO, which is basically a 
standard CDO moved into the future, is that the tranche spreads decrease with the forward start 
date. This is in contrast to the 2007 environment where the risk is generally perceived to increase 
with the forward start date. 

The behaviour is more involved for the non-resetting tranches, where the outstanding tranche 
notional decays leading up to the forward start date. The net effect here is that typically the 
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spreads decrease, with respect to the forward start date, for the junior tranches, increase for the 
senior tranches, and have a hump-shaped term-structure in the mezzanine range. 



6 Conclusions 

We have illustrated the ability of our simple model to crudely calibrate to the index term- 
structure in wildly different market environments, and have shown that the correlation sensitivity 
of tranche spreads demonstrates the behaviour expected. More importantly, using just two param- 
eters and without making them time-dependent, we have shown that our very simple structural 
evolution model displays realistic term-structure dynamics. Using just the volatility parameter, 
it is able to calibrate well to all three index spreads and correlation sensitivities of the various 
tranches are fairly stable across maturities. This is an improvement on the majority of pricing 
models which lack a coherent means of incorporating dynamics. 

The next stage, which has not been the focus here, is to extend the framework so that it can 
calibrate to all tranches with a single set of parameters. This will involve moving away from a 
simple Brownian Motion driving the process, and may include a more general stochastic volatility 
or Levy or jump-diffusion process. Jumps in the market factor are conceptually easy to include 
and result in a jump process driving the SPDE drift. Similarly, a single stochastic volatility factor 
affecting all firms will result in a stochastic term driving the SPDE diffusion. Contagion may 
be incorporated by making model parameters, notably the correlation, loss dependent. These 
extensions allowe the loss distribution process to become more skewed, allocating more weight to 
the tail and increasing super senior tranche spreads, as well as generally allowing more flexibility 
to match observed data. 
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